cpi uiaii 



Integrable impurities as boundary conditions 



M. Moliner 1 ' 2 (a) and P. Schmitteckert 1 ' 2 



1 Institute of Nanotechnology, Karlsruhe Institute of Technology, 76344 Eggenstein-Leopoldshafen, Germany. 

2 Center of Functional Nano structures, Karlsruhe Institute of Technology, 76131 Karlsruhe, Germany. 



PACS 05 . 60 . Gg - Quantum transport 

PACS 73 . 63 . -b - Electronic transport in nanoscale materials and structures 
PACS 73 . 63 . Nm - Quantum wires 

Abstract. - A few exactly solvable interacting quantum many-body problems with impuri- 
ties were previously reported to exhibit unusual features such as non-localization and absence of 
backscattering. In this work we consider the use of these integrable impurities as boundary condi- 
tions in the framework of linear transport problems. We first show that such impurities enhance 
the density of states at the Fermi surface, thus increasing the effective system size. The study 
of the real time-dynamics of a wave packet sent through a series of them inserted in both non- 
interacting and interacting leads then indicates that these impurities are transparent and do not 
add artefacts to the measurement of transport properties. We finally apply these new boundary 
conditions to study the conductance of an interacting scatterer using the embedding method. 



Introduction. — The main difficulty in simulating 
strongly correlated quantum systems consists in the fact 
that strong correlations correspond to low energy, i.e. 
small temperature and large system size which leads to 
exponentially large Hilbert spaces. In his seminal work 
on the solution of the Kondo problem via the Numeri- 
cal Renormalization Group technique (NRG) Wilson [1] 
solved this problem by describing the lead or the bath cou- 
pled to an impurity by an effective model, where he first 
introduced a logarithmic discretization in energy space, 
which he then transformed into a nearest neighbour hop- 
ping chain via a Lanczos tridiagonalization. In this ef- 
fective lead the hopping elements are decaying exponen- 
tially with distance from the impurity leading to an ex- 
ponential increase of the density of states at zero en- 
ergy. This approach has by now been established as 
the reference/benchmark method [2] to solve Kondo like 
problems. Motivated by this success the NRG was re- 
cently extended to describe non-equilibrium properties 
and even steady-state physics of the Single Impurity An- 
derson model (SIAM) [3]. However, in the context of 
transport calculations the Wilson chain possesses one ma- 
jor disadvantage. Each change in the hopping elements 
introduces a scatterer which leads to backscattering, com- 
pare [4]. Thus, one has to be very careful, not to measure 
the transport properties of these "auxiliary" scatterers, 
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especially if one is interested in on-resonance phenomena. 
Here we apply integrable scatterers instead of Wilson's re- 
duced hopping elements which share the property of an 
exponentially enhanced density of states at zero energy, 
while not introducing any backscattering due to their in- 
tegrable nature. 

In order to access larger effective system sizes for bulk 
systems, Vekic and White [5] introduced the Smooth 
Boundary Conditions (SBC) that smoothly reduce the en- 
ergy parameters of the Hamiltonian near the edges of the 
system. Modified boundary conditions were later applied 
to transport problems with the Damped Boundary Con- 
ditions (DBC) [6], where the hopping elements where re- 
duced at the boundaries in a NRG-likc fashion. The idea of 
SBC and DBC is motivated by the logarithmic discretiza- 
tion of the Wilson NRG. However, here one uses the expo- 
nentially decaying hoppings only as a boundary condition 
in order to mimic a larger system size. The DBC are char- 
acterized by two parameters: the damping parameter A, 
t n — t ■ A" , and the length or the region where the hopping 
is reduced, Ld- More precisely, within the L_d/2 sites on 
the right side, the hopping is given by t ■ A™, n = 1 . . . 4p-, 
and similarly on the right side. A restriction of the method 
is the tuning of the parameter A. Indeed, a too strong 
damping amounts to adding another impurity in the sys- 
tem (the DBC region itself) which induces backscattering 
and can give questionable values for the conductance. In 
this work, we propose a new type of boundary conditions 



P-1 



M. Molincr and P. Schmittcckcrt 



Hf = - V sinh 



Uj (2n 3 hj + i — 2n J ^ 1 n : j + i + 2h 3 ^\h. } - 2n.j - (-1)" 3 (c]_ 1 c J - +1 + h.c.)^ 



(V 2 /A + cosh Uj ) (cosh Vj-l) 


-( c k-i 




sinh v 3 




+ i \j 1 ^—{ 2 cosh ^ ( c I+i c i-i 


-h.c.) - 


-v(4 c 3-i- 


L L 
3=1 3=1 




1 1 N " 

3=1 



sinh !A- 



(1) 



(2) 



based on integrable impurities that induces no backscat- 
tering. In the first part, we show that these impurities 
share the property of enhancing the density of states at 
the Fermi surface. The system thus appears to be effec- 
tively larger as long as one studies only properties which 
are dominated by processes at the Fermi surface. Next we 
show that the impurities are transparent for wave packets 
and that this property is preserved in the extension to an 
interacting chain. 

Finally wc apply our boundary condition to study the 
linear transport of an interacting two site impurity via the 
embedding method. [6, 7] 

Integrable impurities. — The integrable impurities 
we are using were previously derived by Schmittcckcrt et 
al. [8,9] using the quantum inverse scattering method [10] 
and the algebraic Bethe ansatz. This method makes it 
possible to determine the eigenvalues of Hamiltonians that 
correspond to specific scattering matrices which fulfil the 
Yang-Baxter equation. These integrable impurities were 
derived in the context of the spin— 1/2 XXZ model which 
can be transformed to a model of spinless fcrmions via a 
Jordan- Wigner transformation [11]. Here we are using the 
simplest type of impurities obtained by this method and 
one impurity consists of three sites {j — 1, j, j + 1} coupled 
with the fine-tuned parameters given in the integrable sin- 
gle impurity Hamiltonian Hj cq. (1) above. V stands for 
the nearest-neighbour repulsion in the spinless fcrmions 
chain in which the impurities are implemented and hj is 
the density operator on site j. The parameter vj char- 
acterizes the strength of the impurity centred on site j. 
Njj integrable impurities are implemented in an interact- 
ing chain of length L with hopping parameter t = 1 as in 
eq. (2) above. 

In non-interacting leads (V = 0) such as the ones used 
in the embedding method, it simplifies as [9] 

L 
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Four of these integrable impurities are represented in the 
boundary region (green) in fig. 3. Setting v 3 = removes 
the impurities and eq. (3) boils down to a one-dimensional 
tight-binding Hamiltonian. On the opposite, Vj — > oo dis- 
connects the site j from the chain and adds a phase e i7r ^ 2 
between the sites j — 1 and j + 1 . These impurities break 
translational invariancc and replace it by another symme- 
try generated by iA,-dcpcndant operators that can be cal- 
culated using higher derivative of the transfer matrices. 
As a consequence, no degeneracies arc lifted, contrary to 
generic impurities. Another noteworthy feature of these 
integrable impurities is that they induce no localization, 
as long as the phase on the (j — 1, j + 1) bonds is ir/2. For 
this precise value, the squared modulus of the wavefunc- 
tion of a system in presence of Nn impurities shows an 
extended eigenstate with reduced amplitude on the central 
sites of the impurities. This behaviour is directly related 
to the transparency of the impurity [8, 12], which is in 
accordance with the findings of [13]. 

Density of states. — We start by looking at the 
change of the density of states (DOS), i.e. the inverse 
level spacing, which is generated in the presence of the 
integrable impurities. We performed exact diagonaliza- 
tions of a non-interacting lead of length L = 150 pierced 
by an Aharonov-Bohm flux $ = 7r/2 to remove degenera- 
cies of the non-current carrying system. In the following, 
we always implement an even number Nu of impurities 
as boundary conditions, half of them with a strength of 
+v and the other half — v so that the ground state en- 
ergy is an even function of the phase (j) induced by the 
Aharonov-Bohm flux. The upper panel of fig. 1 shows the 
influence of the parameter v in a ring containing 12 inte- 
grable impurities. The density of state at the Fermi level 
increases with v, leading to an extremely sharp peak at 
v = ht. In the lower panel, |^| is fixed at the intermediate 
value \i>\ = 2 and we study the influence of the number of 
impurities Njj (i.e. the size if the boundary region). The 
density of states at the Fermi level also increases with Njj. 

Wave packet dynamics. — Wc now consider the 
real time dynamics of a wave packet sent on a series of 
six integrable impurities first in a non-interacting chain of 
spinless fermions eq. (3) and then in an interacting Lut- 
tinger lead eq. (2) with repulsion V — 1.0. The bulk sys- 
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Fig. 1: (Colour on-line) Density of states obtained by exact 
diagonalizations in a L = 150 non-interacting ring pierced by 
a flux $ = 7r/2 with Nn integrable impurities as boundary 
conditions {y — ±H). Upper panel: influence of \v\ at Nu = 
12 fixed. Lower panel: influence of Nn at \v\ = 2.0 fixed. 



Fig. 2: (Colour on-line) Snapshots of the evolution of a Gaus- 
sian wave packet, starting from site 15, through six integrable 
impurities (y — 1) centred on sites j = 31 ... 41 in a non- 
interacting (top panel) and in an interacting system (bottom 
panel). For comparison, the curves with crosses and thin lines 
show the same system at the same time without impurities. 



tern is half- filled (n = 1/2) and the excitation is created 
by a Gaussian wave packet of right-moving fermions with 
wave-vector k = 0.6 (k > kp). We use a time-dependant 
Density Matrix Renormalization Group (td-DMRG) algo- 
rithm based on applying the full time evolution operator 
to an initial state [14], making 11 sweeps and keeping up 
to 1000 states in a system with L = 90 sites and periodic 
boundary conditions (PBC). Fig. 2 shows the evolution of 
the density of the wave packet sent on the series of inte- 
grable impurities with v = 1.0 (thick solid lines) and in 
a system without impurities (curves with crosses) at the 
same time values. Due to band curvature effects the wave 
packet already acquires a side peak in the homogeneous 
system. This feature is preserved in the system including 
the impurities which are indeed completely transparent as 
expected from the Bethe ansatz intcgrability of the model. 
In case of an interacting wire the wave packet splits into 
a dominant right-moving part and a smaller left-moving 
part. The splitting of the bare electron into left and right- 
moving low-lying excitations of the Luttinger liquid is due 
to the density density interaction and the fact, that the 
fermions themselves are not the bare constituents of the 
low-lying excitations. Note that this splitting is also ob- 
served in simulations of the spin charge separation [15] in 
one-dimensional Hubbard chains. Even in the interacting 
case no additional feature appears in the presence of in- 



tegrable impurities, which again illustrates the fact that 
they induce no backscattering. However, since the hop- 
ping parameters are modified in the impurities, the wave 
packet acquires a phase that depends on the value of v 
and on the number of impurities it crosses. In the limit 
v — y oo, the wave-packet would travel trough a smaller 
system of size = L — N n . While this result is to be 
expected from the Bethe ansatz intcgrability of the model 
it displays an example where the intcgrability of the sys- 
tem obtained within an equilibrium description persists in 
the non-equilibrium dynamics. Note that in contrast to 
the NRG tsunami [4] the density of the wave packet is not 
enhanced in IIBC region, see Fig. 2. Thus the IIBC region 
only leads to an increase of the density of states, while the 
modification on the velocity is small. 

The embedding method revisited. Determin- 
ing the conductance of low-dimensional strongly corre- 
lated structures, such as nanotubes or molecules, is now 
an experimental and theoretical challenge [16]. Since it 
is currently not possible to study a realistic model for 
transport through molecules including a complete quan- 
tum chemistry description using a many-particle approach 
we restrict ourselves to the study of model systems of in- 
teracting structures as benchmarks for more complex sys- 
tems. In this section we compute the conductance of a 
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Fig. 3: (Colour on-line) A 4-sites interacting scatterer embed- 
ded in a non-interacting ring pierced by an Aharonov-Bohm 
flux <E>. The boundary conditions (green) consist in integrable 
impurities (here four) with parameter ±v (see eq. (3)). 

small scatterer, namely a few interacting sites, embedded 
in a non-interacting tight-binding ring of length L pierced 
by an Aharonov-Bohm flux <I> (see fig. 3). The embedding 
method [7, 17] is based on the fact that the conductance 
of the scatterer can be related to the persistent current 
generated by the Aharonov-Bohm effect, J = dEo/d$>, 
where Eo is the ground-state energy. In a ring of N spin- 
less fermions (N odd 1 ), with length L and with disorder 
induced by a general potential V, the leading 1/N term 
in the persistent current is [18]: 



Art = - 



evF arccos (|Ifc F cosi 



10 (4) 



where = 2irQ/&o = e&/ch and |T/c F | 2 is the transmis- 
sion probability at the Fermi level. So far, no analytical 
solution was proposed for the case of an interacting scat- 
terer. However, in the limit of large electron reservoirs, 
i.e. leads at least larger than the size of correlations in- 
duced in the lead by the scatterer (e.g. Kondo cloud), and 
under the assumption of Fermi liquid behaviour, the in- 
teracting scatterer behaves effectively as a non-interacting 
scatterer with an interaction-dependent transmission co- 
efficient [19-21]. At very low temperature, the conduc- 
tance of a small interacting region is then obtained by a 
Landauer-like formula [22-24] g = \T kp \ 2 . 

Study of a Kondo-like model. We now use these 
integrable impurities as boundary conditions for the em- 
bedding method to study a Kondo-like model and we com- 
pare the value of the conductance with the one obtained 
by other types of boundary conditions. We consider a non- 
interacting tight-binding ring (bulk region) embedding a 
four-sites interacting region with repulsions U\ and U 2 and 
hopping constants t\ and t 2l as showed in fig. 3. The ring 
is pierced by an Aharonov-Bohm flux $ which is included 
in the first bond and the modified boundary conditions 
are applied between the sites Xi and Xf. 



even adds one additional term to the current but does not 
change the physical content presented hereafter. 
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By looking at the system before closing left and right 
leads to a ring, i.e. in a two-lead set-up, we can label 
the left lead and the first two impurity sites as down spin 
fermions and the other fermions as up spin fermions. By 
merging the up and down sites with the same distance 
from the centre into spin 1/2 sites, c x ,-\ '■= c x and c x ,i '■= 

C— Xl 

we can map the model eq. (5) to an extended single 
impurity Anderson model, see fig. 3. In this description 
U 2 takes the role of the on-site repulsion, t± corresponds 
to the hybridization to the leads, while t 2 plays the role 
of a transverse magnetic field *2(4 ,t c so4 + cl oX c XoA ) ~ 
^{Sq + Sq). t 2 has thus a double role in the trans- 
port properties. Setting t 2 to zero disables any transport, 
while, as detailed later, a finite t 2 reduces and for t 2 > 7k 
even destroys the Kondo peak. Here we are looking at a 
regime, where t 2 is finite to enable transport, yet small 
enough to not destroy the Kondo effect. The repulsion 
Ui(ni ta — 1/2)(77o,,t — 1/2) plays a similar role as in the 
Interacting Resonant Level Model (IRLM) where it broad- 
ens the conductance resonance peaks. [25,26] Hereafter we 
set U 2 = 0.6, t\ = 0.3 and t 2 — 0.05 while the lead hopping 
t is used as an energy scale and U\ varies. We studied the 
effect on the conductance of an asymmetric gate voltage, 
which in the spin language corresponds to a small mag- 
netic field B, to get an approximate value of the Kondo 
temperature which was found to be Tk ~ 0.1. 

In the following we use the DMRG, making 11 sweeps 
and keeping up to 1700 states, to compute the expectation 
value of the current operator Jj = — ij^Im(ctc J -_ 1 ), 

where tj is the hopping matrix element between sites 
j and j — 1. This value is then fitted with eq. (4) to 
extract the conductance g and the effective system size 
L e g = VF/L e ff, where vf is the Fermi velocity. Note that 
L e s should be understood as an effective size only and 
the results may have deviations from currents measured 
in a system of L e fF systems sites in the absence of modi- 
fied boundary terms due to effects of the Kondo screening 
cloud [27]. 

Fig. 4 shows the conductance, extracted from DMRG, 
versus the repulsion U 1 for different types of boundary con- 
ditions: PBC, DBC and Integrable Impurities as Bound- 
ary Conditions (IIBC). DBC (rcsp. IIBC) were applied 
on 24 bonds with Mn = 12 on each side (resp. Nu = 12 
integrable impurities) in L = 30 rings with different values 
of the damping parameter A (rcsp. v = ±|z/|). For com- 
parison, rings with L — 30 and L = 190 sites and PBC 
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Fig. 4: (Colour on-line) Conductance g versus Ui for differ- 
ent types of boundary conditions. DBC (diamonds) and IIBC 
(circles) were applied on systems with L = 30 for different val- 
ues of A and v. For comparison, systems with PBC are also 
shown (crosses). As expected from the definition of the mod- 
ified boundary conditions, one recovers exactly the values of 
the PBC (L = 30) for A = 1.0 and v — (orange curves). 



are also shown. On a qualitative level all three types of 
boundary conditions provide the same result. We first see 
a finite conductance, which is reduced from the unitarity 
limit of g = 2e 2 jh by the scattering processes introduced 
by ti. Interestingly, similar observations were made for the 
IRLM. [25] Interaction on the contact link first increases 
the conductance, even up to the unitarity limit, while for 
large Ui the conductance decreases again. Looking at the 
quantitative results, the different boundary schemes lead 
to different results. To elucidate the differences further we 
compare in fig. 5 for three values of U\ (U\ = 0.0, 0.4, 0.7) 
the conductance versus the effective system size for differ- 
ent types of boundary conditions. PBC were applied in 
systems with up to L = 430 sites (squares), DBC (dia- 
monds) and IIBC (circles) in smaller systems. DBC were 
applied on 12 sites at each end of a L = 30 sites rings with 
varying damping parameters A = 1.0 ... 0.5 (green plots). 
On the one hand, it appears clearly in the centre panel of 
fig. 5 (Ui = 0.4) that for values below A w 0.85, the damp- 
ing tends to reduce the conductance [28], as an effect of 
backscattering on the boundary conditions. For U\ = 0.0 
(resp. U\ = 0.7), decreasing the value of the damping 
parameter A from 1 to 0.5 tends to decrease (resp. in- 
crease) the conductance as the system size increases. On 
the other hand, the size of the damping region at fixed 
A = 0.85 does not have much effect on the conductance 
(curves with purple filled diamonds) [6]. 

IIBC were applied on 24 sites (Njj = 12) in L = 30 sites 
rings with increasing parameter |^| = . . . 10. A sweeping 
procedure, where we gradually switched on the v parame- 
ter, was used in the same way as in [6] to avoid trapping of 
the DMRG in some excited states due to the low coupling 
of the central sites of the impurities for large values of 
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Fig. 5: (Colour on-line) Conductance g versus l/I/ e ff obtained 
from eq. (4) for different types of boundary conditions: PBC 
(squares), DBC (diamonds) and IIBC (circles). Top panel: 
Ui = 0.0. Centre panel: Ui = 0.4. Bottom panel: Ui = 0.7. 



|^| . Fig. 5 shows that increasing |^| first tends to increase 
the effective system size, which is the expected effect, 
and gives values of the conductance comparable with the 
ones obtained with PBC and DBC. However for v > 2.5, 
the conductance and the effective system size strongly de- 
creases, as the central sites of the impurities start to dis- 
connect from the chain. The effective system size finally 
converges to a value related to the size of the ring without 
the impurities central sites L — Nj] , and the conductance 
saturates. We also studied the effect of the number of 
integrable impurities at fixed =2 (orange plots). It 
appears that increasing the size of the boundary region 
tends to increase the conductance, eventually leading to a 
sudden upturn. Actually, even for a single weak link with- 
out any interaction we find that in the large v limit the 
transmission obtained from the embedding method shows 
an upturn toward unitary transmission. We would like to 
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Fig. 6: (Colour on-line) Conductance g versus £2 for a L = 90 
system with Nil = 12 (triangles), 36 (crosses) and 42 (circles) 
integrable impurities with v — |2.0|. The repulsion U\ broadens 
the region of conductance g ~ 1: Ui = 0.0 (orange plots), 
Ui — 0.4 (blue plots) and Ui = 0.7 (green plots). 



note that the embedding method is restricted to strict one- 
dimensional leads in order to avoid the problem of channel 
mixing, for details see [28]. A system with no scatterer is 
integrable via the algebraic Bethe ansatz and we expect 
strict onc-dimcnsional behaviour as given by the analytic 
solution of the ground state curvature [8,9]. However, the 
scatterer under investigation breaks the integrability of 
the system which could finally lead to a channel mixing at 
the impurity and a breakdown of the embedding method. 
We would like to point out, that this is is not a problem 
of our integrable boundary condition; at least not as long 
as one does not extract quantities which rely on a strictly 
one-dimensional character of the system. 

It is important to keep in mind that the conductance 
should be extracted by performing the thermodynamic 
limit to infinite lead length. Ideally, it should be indepen- 
dent from the choice of boundary conditions, as in fig. 5 
top and middle panel. However, as seen in fig. 5 lower 
panel the conductance obtained from different boundaries 
conditions can substantially differ, signalling that the re- 
sult gives only the order magnitude, but not the precise 
quantitative value. By applying different types of bound- 
ary conditions, one has a handle to check the sensitivity of 
a system to these issues. We recommend that one should 
use both types of boundary conditions to check results for 
consistency instead of relying on DBC and IIBC alone. 
Finally, we study the effect of the magnetic field £2 on the 
conductance at the same three values of U\ (fig. 6). As 
expected by the Kondo physics previously discussed, the 
Kondo peak is destroyed by a magnetic field ti <~ Tk- 
Moreover the range of values of t-i for which the conduc- 
tance reaches .9 = 1 becomes wider as \J\ increases (sharp 
peak for XJ\ = 0.0, wide plateau for U\ = 0.7) since this 
interaction tends to broaden the resonance peak. The de- 
crease of the conductance for small ti is related to the 
shrinking of the resonance width while ti goes to zero. 
E.g. for t-i = it has to be strictly zero. Thus, the unitary 
conductance values get reduced by our limited resolution 



stemming from the finite level spacing of the leads. 

Conclusion. — We have proposed a new type of 
boundary conditions based on integrable impurities to in- 
crease the effective system size in the frame of the embed- 
ding method. Unlike other modified boundary conditions, 
such as DBC, these integrable impurities generate no in- 
terfering backscattering. We have first justified their use 
for linear transport problems by showing that they indeed 
increase the effective system size. We have then showed 
that they are transparent for a wave-packet even in an in- 
teracting system. We have applied these new boundary 
conditions to study the transport properties of a system 
that displays geometric Kondo physics and showed that 
the Kondo resonance gets broadened by an interaction on 
the contact link. We would like to emphasize that the use 
of the IIBC is not restricted to DMRG. It is a very general 
concept and can be easily implemented in other methods 
like density functional theories or quantum Monte Carlo 
techniques. By making v site dependent they can also be 
implemented within standard or time dependent NRG. 
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